--------------------------------------------------------------------------------------       log:  /jss/Staiger Productivity/stage1_hosp_nov07.log  log type:  text opened on:  26 Nov 2007, 17:51:43. set mem 800m(819200k). set more off. . * SET GLOBALS -- DEP VAR, CONTROLS, ETC . . global id_hosp "provider". . * 2. the dependent variables -- mortality, part A cost, and total cost . . global depvar "dead1yr". . * 3. the control variables for each cohort: . . global covar_ami "$asrace $riskadj". . /* this global determines the minimum size threshold */. /* threshold = cut prior to 92, and 5*cut 92-04 */. global cut=1. . . set mem 1200m(1228800k). set matsize 800. . use  "/jss/Trend/ami/Unified_ami/unified_ami8604.dta",clear. . drop if mi == 1(141158 observations deleted). drop if i_hmo==1(114236 observations deleted). . * Goes back. *drop if year < 1992. . *  Create an extract. *  keep if 500*int(_n/500) == _n. . tab year, gen(yr)       year |      Freq.     Percent        Cum.------------+-----------------------------------       1986 |     47,113        1.48        1.48       1987 |     45,586        1.43        2.91       1988 |     45,620        1.43        4.34       1989 |     44,401        1.39        5.74       1990 |     44,818        1.41        7.14       1991 |     51,482        1.62        8.76       1992 |    248,927        7.81       16.57       1993 |    240,119        7.54       24.11       1994 |    234,580        7.36       31.47       1995 |    229,568        7.21       38.68       1996 |    224,404        7.04       45.72       1997 |    215,385        6.76       52.48       1998 |    214,328        6.73       59.21       1999 |    216,518        6.80       66.01       2000 |    220,531        6.92       72.93       2001 |    221,489        6.95       79.88       2002 |    224,123        7.03       86.92       2003 |    217,223        6.82       93.73       2004 |    199,622        6.27      100.00------------+-----------------------------------      Total |  3,185,837      100.00. . rename sum_ptareimbrs parta. rename sum_icudays icu. rename sum_ccudays ccu. . . destring dx1, replacedx1 has all characters numeric; replaced as long. gen dxshort = int(dx1/10) if dx1 > 5000(170466 missing values generated). replace dxshort = dx1 if dx1 <= 5000(170466 real changes made). . gen anterlat  = (dxshort==4100). gen antwall   = (dxshort==4101). gen inflat    = (dxshort==4102). gen infpost   = (dxshort==4103). gen infwall   = (dxshort==4104). gen latnec    = (dxshort==4105). gen truepost  = (dxshort==4106). gen subend   = (dxshort==4107). gen aminec    = (dxshort==4108). gen aminos    = (dxshort==4109). . * Create age-sex-race variables. gen iage=recode(age,69,74,79,84,110). egen iasr = group(sex race iage). tab iasr, gen(asr)  group(sex | race iage) |      Freq.     Percent        Cum.------------+-----------------------------------          1 |     24,699        0.78        0.78          2 |     26,454        0.83        1.61          3 |     25,195        0.79        2.40          4 |     20,830        0.65        3.05          5 |     24,868        0.78        3.83          6 |     11,723        0.37        4.20          7 |     11,897        0.37        4.57          8 |     11,490        0.36        4.93          9 |     10,356        0.33        5.26         10 |     11,632        0.37        5.62         11 |    203,541        6.39       12.01         12 |    264,291        8.30       20.31         13 |    303,663        9.53       29.84         14 |    297,957        9.35       39.19         15 |    357,406       11.22       50.41         16 |     24,234        0.76       51.17         17 |     21,622        0.68       51.85         18 |     17,172        0.54       52.39         19 |     11,736        0.37       52.76         20 |      9,910        0.31       53.07         21 |     17,707        0.56       53.62         22 |     14,113        0.44       54.07         23 |     11,567        0.36       54.43         24 |      8,212        0.26       54.69         25 |      6,607        0.21       54.90         26 |    343,819       10.79       65.69         27 |    354,401       11.12       76.81         28 |    316,366        9.93       86.74         29 |    236,890        7.44       94.18         30 |    185,479        5.82      100.00------------+-----------------------------------      Total |  3,185,837      100.00. . *  Collapse a few variables. gen diabetes = (diabet1==1 | diabet2 == 1). gen liver = (liver1==1 | liver2==1). gen cancerall = (cancer==1 | metast==1). . #del;delimiter now ;. global riskadj "vascul pulmon dementia diabetes liver renal cancerall > anterlat antwall inflat infpost infwall latnec truepost subend aminec aminos";. #del cr;delimiter now cr. . . global year04    "yr2-yr19". global asrace "asr2-asr30". global covar_ami "$asrace $riskadj". . * Running regressions to get Hospital-specific effects. . gen dead1yr = (survdays < 365) . gen dead6m  = (survdays < 183). gen dead30d = (survdays <  30). gen dead7d  = (survdays <   7). . * create id variable that is hospital-by-year . egen idx = group(provider year). . compressdxshort was float now intanterlat was float now byteantwall was float now byteinflat was float now byteinfpost was float now byteinfwall was float now bytelatnec was float now bytetruepost was float now bytesubend was float now byteaminec was float now byteaminos was float now byteiage was float now intiasr was float now bytediabetes was float now byteliver was float now bytecancerall was float now bytedead1yr was float now bytedead6m was float now bytedead30d was float now bytedead7d was float now byte. summ $depvar $riskadj age    Variable |       Obs        Mean    Std. Dev.       Min        Max-------------+--------------------------------------------------------     dead1yr |   3185837    .3338966    .4716033          0          1      vascul |   3185837    .0691655    .2537354          0          1      pulmon |   3185837    .1869286    .3898543          0          1    dementia |   3185837    .0295963    .1694709          0          1    diabetes |   3185837    .2460443    .4307047          0          1-------------+--------------------------------------------------------       liver |   3185837    .0026659    .0515631          0          1       renal |   3185837    .0237156    .1521616          0          1   cancerall |   3185837    .0417482    .2000133          0          1    anterlat |   3185837    .0420307    .2006593          0          1     antwall |   3185837    .1650976    .3712687          0          1-------------+--------------------------------------------------------      inflat |   3185837     .029333    .1687381          0          1     infpost |   3185837    .0207537    .1425588          0          1     infwall |   3185837    .1779887    .3825033          0          1      latnec |   3185837    .0212462    .1442041          0          1    truepost |   3185837    .0072188    .0846565          0          1-------------+--------------------------------------------------------      subend |   3185837    .4142403    .4925904          0          1      aminec |   3185837    .0217811    .1459681          0          1      aminos |   3185837    .1003099    .3004128          0          1         age |   3185837    77.07313    7.704767         65        190. tab year, summ($depvar)            |         Summary of dead1yr       year |        Mean   Std. Dev.       Freq.------------+------------------------------------       1986 |   .39746991   .48937988       47113       1987 |   .38643443   .48693743       45586       1988 |   .37553705   .48426657       45620       1989 |    .3629648   .48086023       44401       1990 |   .35041724   .47710594       44818       1991 |   .34953964   .47682921       51482       1992 |   .34173874   .47429345      248927       1993 |   .33468822   .47188234      240119       1994 |   .32832722   .46960558      234580       1995 |   .32403471   .46801407      229568       1996 |   .32311367   .46766676      224404       1997 |    .3270237   .46912709      215385       1998 |   .33322758   .47136822      214328       1999 |   .33763013    .4729028      216518       2000 |   .33229342   .47103663      220531       2001 |   .33415655   .47169583      221489       2002 |   .32892207   .46982265      224123       2003 |   .32843207   .46964397      217223       2004 |   .31954394   .46630108      199622------------+------------------------------------      Total |   .33389656   .47160335     3185837. tab sex        sex |      Freq.     Percent        Cum.------------+-----------------------------------       Fema |  1,606,002       50.41       50.41       Male |  1,579,835       49.59      100.00------------+-----------------------------------      Total |  3,185,837      100.00. tab race       race |      Freq.     Percent        Cum.------------+-----------------------------------      Black |    206,720        6.49        6.49      Other |    115,304        3.62       10.11      White |  2,863,813       89.89      100.00------------+-----------------------------------      Total |  3,185,837      100.00. tab iage       iage |      Freq.     Percent        Cum.------------+-----------------------------------         69 |    625,723       19.64       19.64         74 |    692,778       21.75       41.39         79 |    685,453       21.52       62.90         84 |    585,981       18.39       81.30        110 |    595,902       18.70      100.00------------+-----------------------------------      Total |  3,185,837      100.00. . /* now select the sample */. . sort provider. save temp, replacefile temp.dta saved. . gen nobs=1. summ nobs    Variable |       Obs        Mean    Std. Dev.       Min        Max-------------+--------------------------------------------------------        nobs |   3185837           1           0          1          1. collapse (sum) nobs, by(provider year). gen flag_cut = 0. replace flag_cut = 1 if nobs<$cut & year<1992(0 real changes made). replace flag_cut = 1 if nobs<(5*$cut) & year>=1992(6790 real changes made). egen min_nobs_cut = max(flag_cut), by(provider). gen one=1. egen tot_yrs = sum(one), by(provider). gen insamp = (min_nobs_cut==0 & tot_yrs==19). tab year insamp           |        insamp      year |         0          1 |     Total-----------+----------------------+----------      1986 |     2,306      2,819 |     5,125       1987 |     2,211      2,819 |     5,030       1988 |     2,102      2,819 |     4,921       1989 |     1,978      2,819 |     4,797       1990 |     1,881      2,819 |     4,700       1991 |     1,818      2,819 |     4,637       1992 |     2,257      2,819 |     5,076       1993 |     2,166      2,819 |     4,985       1994 |     2,129      2,819 |     4,948       1995 |     2,065      2,819 |     4,884       1996 |     1,997      2,819 |     4,816       1997 |     1,896      2,819 |     4,715       1998 |     1,860      2,819 |     4,679       1999 |     1,794      2,819 |     4,613       2000 |     1,639      2,819 |     4,458       2001 |     1,429      2,819 |     4,248       2002 |     1,193      2,819 |     4,012       2003 |     1,023      2,819 |     3,842       2004 |       907      2,819 |     3,726 -----------+----------------------+----------     Total |    34,651     53,561 |    88,212 . tab year insamp [fw=nobs]           |        insamp      year |         0          1 |     Total-----------+----------------------+----------      1986 |    10,647     36,466 |    47,113       1987 |     9,731     35,855 |    45,586       1988 |     9,027     36,593 |    45,620       1989 |     8,377     36,024 |    44,401       1990 |     8,197     36,621 |    44,818       1991 |     8,618     42,864 |    51,482       1992 |    40,847    208,080 |   248,927       1993 |    37,725    202,394 |   240,119       1994 |    34,879    199,701 |   234,580       1995 |    32,267    197,301 |   229,568       1996 |    28,789    195,615 |   224,404       1997 |    24,920    190,465 |   215,385       1998 |    22,944    191,384 |   214,328       1999 |    21,646    194,872 |   216,518       2000 |    19,462    201,069 |   220,531       2001 |    17,478    204,011 |   221,489       2002 |    15,490    208,633 |   224,123       2003 |    14,336    202,887 |   217,223       2004 |    12,286    187,336 |   199,622 -----------+----------------------+----------     Total |   377,666  2,808,171 | 3,185,837 . keep if insamp==1(34651 observations deleted). tab year       year |      Freq.     Percent        Cum.------------+-----------------------------------       1986 |      2,819        5.26        5.26       1987 |      2,819        5.26       10.53       1988 |      2,819        5.26       15.79       1989 |      2,819        5.26       21.05       1990 |      2,819        5.26       26.32       1991 |      2,819        5.26       31.58       1992 |      2,819        5.26       36.84       1993 |      2,819        5.26       42.11       1994 |      2,819        5.26       47.37       1995 |      2,819        5.26       52.63       1996 |      2,819        5.26       57.89       1997 |      2,819        5.26       63.16       1998 |      2,819        5.26       68.42       1999 |      2,819        5.26       73.68       2000 |      2,819        5.26       78.95       2001 |      2,819        5.26       84.21       2002 |      2,819        5.26       89.47       2003 |      2,819        5.26       94.74       2004 |      2,819        5.26      100.00------------+-----------------------------------      Total |     53,561      100.00. tab year insamp [fw=nobs]           |   insamp      year |         1 |     Total-----------+-----------+----------      1986 |    36,466 |    36,466       1987 |    35,855 |    35,855       1988 |    36,593 |    36,593       1989 |    36,024 |    36,024       1990 |    36,621 |    36,621       1991 |    42,864 |    42,864       1992 |   208,080 |   208,080       1993 |   202,394 |   202,394       1994 |   199,701 |   199,701       1995 |   197,301 |   197,301       1996 |   195,615 |   195,615       1997 |   190,465 |   190,465       1998 |   191,384 |   191,384       1999 |   194,872 |   194,872       2000 |   201,069 |   201,069       2001 |   204,011 |   204,011       2002 |   208,633 |   208,633       2003 |   202,887 |   202,887       2004 |   187,336 |   187,336 -----------+-----------+----------     Total | 2,808,171 | 2,808,171 . keep if year==1994(50742 observations deleted). keep provider. sort provider. save id_insamp, replacefile id_insamp.dta saved. . . use temp, clear. merge provider using id_insampvariable provider does not uniquely identify observations in the master data. tab _merge     _merge |      Freq.     Percent        Cum.------------+-----------------------------------          1 |    377,666       11.85       11.85          3 |  2,808,171       88.15      100.00------------+-----------------------------------      Total |  3,185,837      100.00. keep if _merge==3(377666 observations deleted). drop _merge. . .  .       /* 1. RUN THE BASIC REGRESSION  */.       .       areg $depvar $covar_ami, absorb(idx)Linear regression, absorbing indicators                Number of obs = 2808171                                                       F( 45,2754565) = 6937.78                                                       Prob > F      =  0.0000                                                       R-squared     =  0.1299                                                       Adj R-squared =  0.1130                                                       Root MSE      =  .44282------------------------------------------------------------------------------     dead1yr |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]-------------+----------------------------------------------------------------        asr2 |    .046187   .0042677    10.82   0.000     .0378225    .0545515        asr3 |   .1021323   .0043212    23.63   0.000     .0936628    .1106018        asr4 |   .1584622   .0045532    34.80   0.000      .149538    .1673863        asr5 |   .2737937   .0043658    62.71   0.000     .2652369    .2823504        asr6 |  -.0283706   .0054582    -5.20   0.000    -.0390685   -.0176727        asr7 |   .0163674   .0054377     3.01   0.003     .0057098    .0270251        asr8 |   .0674053   .0054979    12.26   0.000     .0566296    .0781809        asr9 |   .1582632    .005699    27.77   0.000     .1470933    .1694331       asr10 |    .267205   .0054661    48.88   0.000     .2564916    .2779184       asr11 |  -.0334198   .0032794   -10.19   0.000    -.0398474   -.0269922       asr12 |    .013856   .0032408     4.28   0.000     .0075041    .0202079       asr13 |   .0724861   .0032256    22.47   0.000      .066164    .0788083       asr14 |   .1492104   .0032307    46.19   0.000     .1428784    .1555424       asr15 |   .2790713   .0032176    86.73   0.000     .2727649    .2853776       asr16 |  -.0215054   .0043575    -4.94   0.000     -.030046   -.0129649       asr17 |   .0303059   .0044874     6.75   0.000     .0215106    .0391011       asr18 |   .0897595    .004795    18.72   0.000     .0803614    .0991575       asr19 |    .165311   .0054435    30.37   0.000     .1546419    .1759802       asr20 |   .2775689   .0058094    47.78   0.000     .2661827    .2889551       asr21 |  -.0523235   .0047799   -10.95   0.000     -.061692    -.042955       asr22 |   .0046436   .0051501     0.90   0.367    -.0054504    .0147375       asr23 |   .0748653   .0055101    13.59   0.000     .0640658    .0856649       asr24 |   .1564858   .0062347    25.10   0.000     .1442659    .1687057       asr25 |   .2930539   .0067744    43.26   0.000     .2797764    .3063314       asr26 |  -.0580441   .0032113   -18.07   0.000    -.0643382     -.05175       asr27 |  -.0031456    .003209    -0.98   0.327     -.009435    .0031438       asr28 |   .0713866   .0032222    22.15   0.000     .0650712    .0777021       asr29 |   .1615914   .0032618    49.54   0.000     .1551984    .1679844       asr30 |   .2965322   .0033091    89.61   0.000     .2900465    .3030178      vascul |   .0281913   .0010507    26.83   0.000      .026132    .0302506      pulmon |    .080567   .0006942   116.06   0.000     .0792064    .0819275    dementia |   .1347246   .0016022    84.09   0.000     .1315843    .1378649    diabetes |   .0391664    .000631    62.07   0.000     .0379297    .0404031       liver |   .2399882    .005122    46.85   0.000     .2299492    .2500272       renal |   .2780818    .001784   155.87   0.000     .2745852    .2815783   cancerall |   .1636042   .0013325   122.78   0.000     .1609926    .1662157    anterlat |   .0028791   .0023276     1.24   0.216     -.001683    .0074411     antwall |  -.0325738   .0020276   -16.07   0.000    -.0365477   -.0285998      inflat |  -.0620783   .0024801   -25.03   0.000    -.0669392   -.0572174     infpost |  -.0734464   .0026571   -27.64   0.000    -.0786543   -.0682386     infwall |  -.1017416   .0020207   -50.35   0.000    -.1057022   -.0977811      latnec |  -.0659159   .0026563   -24.82   0.000    -.0711222   -.0607097    truepost |   -.073262   .0036923   -19.84   0.000    -.0804988   -.0660252      subend |  -.1194683   .0019629   -60.86   0.000    -.1233154   -.1156211      aminec |  (dropped)      aminos |   .0931605   .0021291    43.75   0.000     .0889874    .0973335       _cons |    .266698   .0036524    73.02   0.000     .2595395    .2738565-------------+----------------------------------------------------------------         idx | F(53560, 2754565) =      1.518   0.000       (53561 categories).       .       /* save the root mean square error of the regression */.       .        gen rmse_hosp_ami = e(rmse) .        .       /* save the FE, which is the risk-adjusted estimate */.       .        predict fe_hosp_ami if e(sample), d .        .       /* saves the # obs on which FE was based */.       .        egen nobs_hosp_ami = count(fe_hosp_ami) if fe_hosp_ami !=., by(idx) .        . * COSTS. * First, adjust prices using a very simple neanderthal approach. . gen     adj =  .712 if year == 1986(2771705 missing values generated). replace adj =  .732 if year == 1987(35855 real changes made). replace adj =  .757 if year == 1988(36593 real changes made). replace adj =  .786 if year == 1989(36024 real changes made). replace adj =  .816 if year == 1990(36621 real changes made). replace adj =  .844 if year == 1991(42864 real changes made). replace adj =  .864 if year == 1992(208080 real changes made). replace adj =  .884 if year == 1993(202394 real changes made). replace adj =  .903 if year == 1994(199701 real changes made). replace adj =  .921 if year == 1995(197301 real changes made). replace adj =  .939 if year == 1996(195615 real changes made). replace adj =  .954 if year == 1997(190465 real changes made). replace adj =  .965 if year == 1998(191384 real changes made). replace adj =  .979 if year == 1999(194872 real changes made). replace adj = 1.000 if year == 2000(201069 real changes made). replace adj = 1.024 if year == 2001(204011 real changes made). replace adj = 1.041 if year == 2002(208633 real changes made). replace adj = 1.060 if year == 2003(202887 real changes made). replace adj = 1.082 if year == 2004(187336 real changes made). . * Norm to year 2004. replace adj = adj/1.082(2808171 real changes made). replace parta1y  = parta1y/adjparta1y was long now double(2588709 real changes made). .  /* 2.  RUN THE COST REGRESSION  */.       .       areg parta1y $covar_ami, absorb(idx)Linear regression, absorbing indicators                Number of obs = 2808170                                                       F( 45,2754564) = 1443.39                                                       Prob > F      =  0.0000                                                       R-squared     =  0.1140                                                       Adj R-squared =  0.0968                                                       Root MSE      =   23288------------------------------------------------------------------------------     parta1y |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]-------------+----------------------------------------------------------------        asr2 |  -151.6695    224.435    -0.68   0.499    -591.5543    288.2153        asr3 |  -1818.415   227.2515    -8.00   0.000     -2263.82    -1373.01        asr4 |  -4540.187   239.4509   -18.96   0.000    -5009.502   -4070.872        asr5 |   -8150.57   229.5944   -35.50   0.000    -8600.567   -7700.573        asr6 |  -504.5058   287.0445    -1.76   0.079    -1067.103    58.09127        asr7 |   80.85709   285.9644     0.28   0.777    -479.6232    641.3373        asr8 |  -1354.912    289.131    -4.69   0.000    -1921.599   -788.2254        asr9 |   -4645.39   299.7093   -15.50   0.000     -5232.81    -4057.97       asr10 |  -8925.129    287.461   -31.05   0.000    -9488.543   -8361.716       asr11 |  -1011.489   172.4643    -5.86   0.000    -1349.513   -673.4646       asr12 |  -852.0507   170.4327    -5.00   0.000    -1186.093   -518.0086       asr13 |  -1611.165   169.6347    -9.50   0.000    -1943.643   -1278.687       asr14 |  -4521.182   169.8997   -26.61   0.000     -4854.18   -4188.185       asr15 |  -8691.832    169.211   -51.37   0.000     -9023.48   -8360.184       asr16 |  -2104.132   229.1589    -9.18   0.000    -2553.276   -1654.989       asr17 |  -1618.842   235.9926    -6.86   0.000     -2081.38   -1156.305       asr18 |  -2134.084   252.1666    -8.46   0.000    -2628.321   -1639.846       asr19 |   -4078.61   286.2729   -14.25   0.000    -4639.695   -3517.525       asr20 |  -7067.139   305.5125   -23.13   0.000    -7665.932   -6468.345       asr21 |  -1370.252   251.3742    -5.45   0.000    -1862.937   -877.5677       asr22 |   22.19019    270.839     0.08   0.935    -508.6447     553.025       asr23 |  -117.8623   289.7719    -0.41   0.684     -685.805    450.0805       asr24 |  -3485.034   327.8822   -10.63   0.000    -4127.672   -2842.397       asr25 |  -7972.919   356.2603   -22.38   0.000    -8671.176   -7274.661       asr26 |  -1747.958   168.8827   -10.35   0.000    -2078.962   -1416.954       asr27 |  -469.0378   168.7571    -2.78   0.005    -799.7959   -138.2797       asr28 |   -767.995   169.4541    -4.53   0.000    -1100.119   -435.8709       asr29 |  -3434.365   171.5359   -20.02   0.000     -3770.57   -3098.161       asr30 |   -7793.56   174.0224   -44.78   0.000    -8134.638   -7452.483      vascul |   1658.023   55.25467    30.01   0.000     1549.726     1766.32      pulmon |   1367.545   36.50685    37.46   0.000     1295.993    1439.097    dementia |   -4919.82   84.26058   -58.39   0.000    -5084.968   -4754.672    diabetes |   1948.031   33.18281    58.71   0.000     1882.994    2013.068       liver |   -2019.96   269.3647    -7.50   0.000    -2547.905   -1492.014       renal |   1255.941   93.81998    13.39   0.000     1072.057    1439.825   cancerall |  -2715.698   70.07327   -38.76   0.000    -2853.039   -2578.357    anterlat |   684.8354   122.4091     5.59   0.000     444.9179    924.7529     antwall |    1117.44   106.6289    10.48   0.000     908.4507    1326.428      inflat |   331.2845   130.4269     2.54   0.011     75.65229    586.9166     infpost |   1032.925   139.7372     7.39   0.000     759.0453    1306.805     infwall |   332.5648   106.2686     3.13   0.002     124.2822    540.8475      latnec |  -114.2438   139.6921    -0.82   0.413    -388.0353    159.5478    truepost |   825.4786   194.1759     4.25   0.000     444.9007    1206.056      subend |   2091.509    103.226    20.26   0.000      1889.19    2293.828      aminec |  (dropped)      aminos |  -1782.997   111.9705   -15.92   0.000    -2002.456   -1563.539       _cons |   24694.68   192.0754   128.57   0.000     24318.22    25071.14-------------+----------------------------------------------------------------         idx | F(53560, 2754564) =      5.137   0.000       (53561 categories).       .       /* save the root mean square error of the regression */.       .        gen rmse_cost = e(rmse) .        .       /* save the FE, which is the risk-adjusted estimate */.       .        predict hosp_cost if e(sample), d (1 missing value generated). . /* 2.  RUN THE DRG REGRESSION  */.       .       areg drgwgt1y $covar_ami, absorb(idx)Linear regression, absorbing indicators                Number of obs = 2808162                                                       F( 45,2754556) = 2125.30                                                       Prob > F      =  0.0000                                                       R-squared     =  0.0795                                                       Adj R-squared =  0.0616                                                       Root MSE      =  3.3797------------------------------------------------------------------------------    drgwgt1y |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]-------------+----------------------------------------------------------------        asr2 |  -.1056915   .0325711    -3.24   0.001    -.1695297   -.0418533        asr3 |  -.3507815   .0329798   -10.64   0.000    -.4154209   -.2861422        asr4 |  -.7652046   .0347503   -22.02   0.000    -.8333139   -.6970953        asr5 |  -1.292053   .0333203   -38.78   0.000     -1.35736   -1.226746        asr6 |   .0167569   .0416573     0.40   0.687    -.0648899    .0984037        asr7 |   .0330809   .0415005     0.80   0.425    -.0482587    .1144205        asr8 |  -.2382828   .0419601    -5.68   0.000    -.3205231   -.1560425        asr9 |  -.7182024   .0434953   -16.51   0.000    -.8034516   -.6329532       asr10 |   -1.35475   .0417177   -32.47   0.000    -1.436516   -1.272985       asr11 |    .026322   .0250289     1.05   0.293    -.0227337    .0753777       asr12 |  -.0202516    .024734    -0.82   0.413    -.0687294    .0282263       asr13 |  -.2094345   .0246182    -8.51   0.000    -.2576853   -.1611837       asr14 |  -.6938467   .0246567   -28.14   0.000    -.7421729   -.6455205       asr15 |  -1.371988   .0245567   -55.87   0.000    -1.420118   -1.323858       asr16 |  -.2228585   .0332566    -6.70   0.000    -.2880404   -.1576767       asr17 |  -.2340243   .0342484    -6.83   0.000      -.30115   -.1668987       asr18 |  -.3663492   .0365956   -10.01   0.000    -.4380753    -.294623       asr19 |  -.6760372   .0415453   -16.27   0.000    -.7574645   -.5946098       asr20 |   -1.10331   .0443374   -24.88   0.000    -1.190209    -1.01641       asr21 |    .040619   .0364806     1.11   0.266    -.0308818    .1121198       asr22 |   .1040208   .0393055     2.65   0.008     .0269835    .1810582       asr23 |  -.0266352   .0420531    -0.63   0.526    -.1090578    .0557874       asr24 |  -.5326732   .0475838   -11.19   0.000    -.6259358   -.4394105       asr25 |  -1.188803   .0517022   -22.99   0.000    -1.290137   -1.087468       asr26 |   .0358674   .0245091     1.46   0.143    -.0121695    .0839043       asr27 |   .1007752   .0244909     4.11   0.000      .052774    .1487765       asr28 |  -.0347381    .024592    -1.41   0.158    -.0829376    .0134614       asr29 |  -.5059131   .0248941   -20.32   0.000    -.5547047   -.4571215       asr30 |  -1.213705    .025255   -48.06   0.000    -1.263204   -1.164206      vascul |   .2992943   .0080189    37.32   0.000     .2835777     .315011      pulmon |   .2073888   .0052981    39.14   0.000     .1970048    .2177729    dementia |  -.8537341   .0122283   -69.82   0.000    -.8777012   -.8297671    diabetes |   .3459231   .0048157    71.83   0.000     .3364846    .3553617       liver |  -.4322708   .0390915   -11.06   0.000    -.5088888   -.3556528       renal |    .050855   .0136156     3.74   0.000     .0241689    .0775412   cancerall |  -.4017425   .0101694   -39.51   0.000    -.4216741   -.3818109    anterlat |   .0745366   .0177646     4.20   0.000     .0397186    .1093546     antwall |   .1658781   .0154745    10.72   0.000     .1355486    .1962075      inflat |   .0482342   .0189282     2.55   0.011     .0111355    .0853329     infpost |   .1287449   .0202793     6.35   0.000     .0889981    .1684917     infwall |   .0642399   .0154222     4.17   0.000      .034013    .0944669      latnec |  -.0308287   .0202728    -1.52   0.128    -.0705627    .0089052    truepost |   .1186174   .0281797     4.21   0.000     .0633861    .1738487      subend |   .3493444   .0149807    23.32   0.000     .3199828     .378706      aminec |  (dropped)      aminos |  -.3323765   .0162497   -20.45   0.000    -.3642254   -.3005277       _cons |   4.419387   .0278749   158.54   0.000     4.364753    4.474021-------------+----------------------------------------------------------------         idx | F(53560, 2754556) =      2.581   0.000       (53561 categories).       .       gen rmse_drg = e(rmse).       .        predict hosp_drg if e(sample), d (9 missing values generated). . sum $riskadj if year <=2004    Variable |       Obs        Mean    Std. Dev.       Min        Max-------------+--------------------------------------------------------      vascul |   2808171    .0704722    .2559412          0          1      pulmon |   2808171    .1887435    .3913048          0          1    dementia |   2808171    .0294053    .1689396          0          1    diabetes |   2808171    .2468511    .4311794          0          1       liver |   2808171     .002716    .0520445          0          1-------------+--------------------------------------------------------       renal |   2808171    .0232407    .1506673          0          1   cancerall |   2808171    .0425693    .2018842          0          1    anterlat |   2808171    .0409113    .1980848          0          1     antwall |   2808171    .1633772    .3697095          0          1      inflat |   2808171    .0288782    .1674643          0          1-------------+--------------------------------------------------------     infpost |   2808171    .0212224    .1441249          0          1     infwall |   2808171    .1775408    .3821258          0          1      latnec |   2808171    .0211273    .1438086          0          1    truepost |   2808171    .0071947     .084516          0          1      subend |   2808171    .4258704    .4944744          0          1-------------+--------------------------------------------------------      aminec |   2808171    .0209952     .143368          0          1      aminos |   2808171    .0928825    .2902678          0          1. .      .    /* collapse and save the data */.   collapse (mean) dead1yr rmse_* hosp_cost parta1y hosp_drg drgwgt1y fe_* nobs_*, by> (provider year).   sort provider.   save fe_ami_hosp_86_nov07, replace(note: file fe_ami_hosp_86_nov07.dta not found)file fe_ami_hosp_86_nov07.dta saved. . set more on. log close       log:  /jss/Staiger Productivity/stage1_hosp_nov07.log  log type:  text closed on:  26 Nov 2007, 18:27:16--------------------------------------------------------------------------------------